Galactic Substructure and Energetic Neutrinos from the Sun and the Earth 
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We consider the effects of Galactic substructure on energetic neutrinos from annihilation of weakly- 
interacting massive particles (WIMPs) that have been captured by the Sun and Earth. Substructure 
gives rise to a time- varying capture rate and thus to time variation in the annihilation rate and result- 
ing energetic-neutrino flux. However, there may be a time lag between the capture and annihilation 
rates. The energetic-neutrino flux may then be determined by the density of dark matter in the 
Solar System's past trajectory, rather than the local density. The signature of such an effect may 
be sought in the ratio of the direct- to indirect-detection rates. 
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Numerous experimental probes have confirmed in- 
directly the presence of a yet unknown form of 
gravitationally-interacting matter in Galactic halos that 
contributes roughly 20% of the total cosmic energy den- 
sity. It is generally assumed that "dark matter" is in 
the form of some yet undiscovered elementary parti- 
cle. Among the plethora of proposed theoretical par- 
ticle dark-matter candidates, weakly interacting mas- 
sive particles (WIMPs) are favored because they provide, 
quite generally, the correct relic abundance and because 
they may be experimentally accessible in the near fu- 
ture. WIMPs arise naturally in supersymmetric exten- 
sions (SUSY) of the Standard Model [T] as well as in 
models with Universal Extra Dimensions (UEDs) [2J. 

The two principle avenues toward dark-matter detec- 
tion are direct detection (DD) via observation of the re- 
coil of a nucleus, when struck by a halo WIMP, in a low- 
background experiment [3| H] ; and neutrino indirect de- 
tection (idD) via observation of energetic neutrinos from 
annihilation of WIMPs that have been captured in the 
Sun (and/or Earth) [SJE]. 

The DD rate is proportional to the local dark-matter 
density. The v\T> rate is proportional to the rate at which 
WIMPs annihilate in the Sun, which in turn depends on 
an integral of the square of the dark-matter density over 
the volume of the Sun. However, the WIMPs depleted in 
the Sun by annihilation are replenished by the capture 
of new WIMPs. In most cases where the vlD signal is 
large enough to be detectable, the timescale for equilibra- 
tion of capture and annihilation is small compared with 
the age of the Solar System. The i^ID rate is then also 
determined by the local dark-matter density. Since the 
capture rate is controlled by the same elastic-scattering 
process that occurs in DD, the DD and i/ID rates are 
roughly proportional [7]. 

In this Letter we investigate the effects of Galactic sub- 
structure on this canonical scenario. Analytic arguments 
and numerical simulations suggest that realistic Galactic 
halos should have significant substructure, remnants of 
smaller halos produced in early stages of the structure- 
formation hierarchy (which may themselves house rem- 



nants of even smaller structures, and so on) [8]. Theoreti- 
cal arguments suggest that the substructure may be scale 
invariant [S] with subhalos extending all the way down to 
sub-Earth-mass scales [TU]. The local dark-matter den- 
sity of different locations at similar Galactocentric radii 
in the Milky Way may thus differ by a few orders of 
magnitude. The analytic descriptions of substructure are 
rough, and the simulations are limited by finite resolu- 
tion, and this motivates the pursuit of avenues toward 
empirically probing the existence of substructure. 

The purpose of this Letter to show that measurements 
of the ratio of DD to idD rates can be used to test for 
Galactic substructure. If there is Galactic substructure, 
then the dark-matter density at the position of the So- 
lar System may vary with time. There is a finite time 
lag between capture and annihilation, and so the cur- 
rent energetic-neutrino flux may be determined not by 
the local dark-matter density, but rather by the density 
of dark matter along the past trajectory of the Solar Sys- 
tem. The ratio for the idD/DD rate may thus differ from 
the canonical prediction. A departure from the canonical 
ratio would thus, if observed, provide information about 
Galactic substructure. Since the equilibration timescale 
in the Earth is generally different from that in the Sun, 
additional information might be provided by observation 
of energetic neutrinos from WIMP annihilation in the 
Earth. 

To illustrate, we suppose the WIMP has a scalar cou- 
pling to nuclei, but the formalism can be easily general- 
ized to spin-dependent WIMPs. Then the DD rate for a 
WIMP of mass m x from a target nucleus of mass m, is 

i% c D = 2.2 x 10 s kg- 1 yr- 1 ft e ,o.»'?eK l m i ) ( i <^ e v ) 

x (locT^ev) (^~T^) a40 ' (1 

where p x ,o.3 is the local dark-matter density in units 
of 0.3 GeV cm~ 3 , and rj c (m x ,mi) (given in Ref. [4]) 
accounts for form-factor suppression. Here, a 40 is the 
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cross section for WIMP-nucleon scattering in units of 
lCT 40 cm 2 . 

The flux of upward muons induced in a neutrino tele- 
scope by neutrinos from WIMP annihilation in the Sun 
is 

r„,„ - 7.3 x 10 5 km~ 2 yr- 1 (iV/iV cq ) 2 p x , .3 

x /(m x )K(m x )/0.1](m x /100GeV) 2 a 40 , (2) 



while the corresponding flux from the Earth is obtained 
by replacing the prefactor of Eq. ^ by 15km -2 yr -1 . 
The function f(m x ) varies over the range 5 < /(m x ) < 
0.5 over the mass range 10 < m x /GeV < 1000 for the 
Sun (with a slightly larger range for the Earth), while 
the function £(to x ) is in the range ~ 0.01 — 0.3 over the 
same mass range. 

The factor N/N cq in Eq. ^ quantifies the number of 
WIMPs in the Sun. Once WIMPs are captured in the 
Sun, they accumulate deep within the solar core, where 
they may annihilate to a variety of heavy Standard Model 
particles which then decay to produce high-energy neu- 
trinos (which may escape the Sun) . The time t evolution 
of the number N of WIMPs in the Sun, is governed by 
the differential equation, 



dN/dt = C c - C a 7V 2 



(3) 



where C c is the capture rate of WIMPs by the Sun, and 
C a N 2 is twice (because each annihilation destroys two 
WIMPs) the effective annihilation rate. If both C c and 
C a are constant and the initial condition is N(t = 0) = 
iVo, the solution to this equation is 



N(t) 



7e" 



-t/r 



C a e*/r + 7 e-*/r ' 



where 
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< 1. 
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and t = l/\/C c C a is the equilibration timecale. After a 
time t>T, the number N approaches N cq = N(t ^S> t) = 
yj C c /C a , and the annihilation rate T v becomes equal to 
(one half) the capture rate, T u = C a N 2 /2 = C c /2. 

The capture rate C c and annihilation coefficient C a , 
and thus the equilibration timescale r, are determined by 
the cross sections for WIMPs to annihilate and to scatter 
from nuclei. The equilibration timescale evaluates to 



r© = 1.9 x 10 5 yrs[p Xi0 . 3 /(m x )(a A w) 26 ] 1 ' 2 
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Here, (04^)26 is the annihilation cross section (times 
relative velocity v in the limit v — > 0) in units of 
1Q -26 cm _3 s _1 . The equilibration timescale for the 
Earth is obtained by replacing the prefactor of Eq. ([6| 
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FIG. 1: The flux of energetic neutrinos from the Sun versus 
the rate for direct detection. Each point denotes a supersym- 
metric model with the correct relic density and consistent with 
experiment. The different symbols indicate the timescale for 
equilibration between capture and annihilation in the Sun. 
The horizontal line indicates a flux-threshold target for fu- 
ture vID experiments (IceCube+DeepCore [13] and the verti- 
cal line a rate-threshold target for DD in a 3-ton liquid-xenon 
detector [14]. For current bounds see 1171. 



by 1.1 x 10 8 yr. Using the canonical numbers we have 
adopted, the equilibration timescales for the Sun and 
Earth are both small compared with the age of the So- 
lar System, but the equilibration timescales may vary 
by several orders of magnitude over reasonable ranges 
of the WIMP parameter space (and even more if more 
exotic physics, like a Sommerfcld |11] or self-capture en- 
hancement [12] . is introduced). To illustrate, we show in 
Fig. [l] the equilibration timescales, for various DD and 
idD rates, for realistic supersymmetric dark-matter can- 
didates (using DarkSUSY [TS]). 

Suppose now the WIMP model parameters are deter- 
mined, e.g., from the LHC and/or by theoretical assump- 
tion/modeling. Then the unknowns in Eqs. (JT|) and ^ 
will be the halo density p x and the number N of WIMPs 
in the Sun (or Earth). The measured DD rate will then 
provide the local halo density p x . Measurement of the 
^ID rate will then determine (N/N cq ) = tanh(t/r) (in 
both the Sun and the Earth). 

For example, suppose the equilibration timescale is 
t w 10 7 years in the Sun, and that the Solar Sys- 
tem entered a region of density p = 100 pq (where 
Pq = 0.3 GeV cm~ 3 is the smoothed local halo den- 
sity) a time St w 10 6 years ago, e.g., a 1O 9 M0 halo 
(see Fig. [2j. We would then see a boosted DD rate 
and a boosted energetic-neutrino flux from the Sun, but 
the energetic-neutrino flux from the Earth would be 
correspondingly weaker, since the equilibration time in 
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FIG. 2: The neutrino-flux enhancement from an encounter, 
of duration 10 7 yr, of the Solar system with a region where 
the dark-matter density is enhanced by a factor of 100 (e.g. 
a 10 9 Mq subhalo). The top panel shows the capture rate 
(i.e., the DD rate). The bottom panel shows the resulting 
energetic-neutrino flux from such an encounter for three equi- 
libration timescales. 



FIG. 3: The neutrino flux in a hypothetical scenario where 
all dark matter is in dense objects of 1M©. Different curves 
correspond to equilibration timescales as shown. Short equi- 
libration timescales (e.g. Sun) almost deplete completely the 
amount of WIMPs in a time comparable to the timescale 
between interactions. Longer equilibration timescales (e.g. 
Earth) result in a constant elevated flux. 



the Earth is longer. Now, suppose that the Solar Sys- 
tem exited this high-density region a million years ago. 
The DD rate would be at the canonical value, but the 
energetic-neutrino fluxes from the Sun/Earth would still 
be boosted. Finally, suppose that the Solar System ex- 
ited the high-density region 10 7 years ago. In that case, 
the DD rate and energetic-neutrino flux from the Sun 
would have the canonical values, but the neutrino flux 
from the Earth would still be boosted, as r e > t . 

In reality, the capture rate C c (t) in Eq. ([3| is a function 
of time, and the equation for the number of WIMPs in 
the Sun or Earth can be integrated numerically to give 
the annihilation rate C a N 2 (t) /2 as a function of time. 
To illustrate, imagine that all of halo dark matter was 
distributed in objects of a single mass, M — M\Mq, 
each with a density 1000 /3p Q [16]. The radius of these 
subhalos would then be R = 0.1 pc (M 1 /(3) 1 / 3 . The tran- 
sit time of the Solar System through such an object is 
St « 1000(Mi//3) 1 / 3 yr. The mean-free time between en- 
counters with such objects is t « 10 7 (Mi/3 2 ) 1 / 3 yr. In 
this toy model, the dark-matter density (and hence the 
DD rate) is zero unless the Solar System is within a sub- 
halo. Fig. [3] illustrates the effects from such a scenario 
for Mi = 1. For equilibration timescales which are of 
order t ~ t , the energetic-neutrino signal is depleted 
completely prior to the next encounter, while for longer 
equilibration timescales, the net effect is an elevated sig- 
nal at all times. For example, if r Q ~ t, and r ffi > r Q , 
the signal from the Earth will be boosted relative to the 
signal from the Sun for most of the time. 
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FIG. 4: The net effect on the energetic-neutrino flux of the 
presence of 1O _6 M0 objects along the solar Galactic radius 
during the whole lifetime of the Solar System. Long equili- 
bration timescales result in the net build-up of WIMPs and 
thus an increase in the neutrino flux relative to that which 
would be obtained in a smooth halo. 

Finally, substructure may also speed up the equilibra- 
tion between capture and annihilation in cases where the 
smooth-halo equilibration timescale is larger than the age 
of the Solar System. Suppose the dark matter has a 
smooth component and some substructure down to very 
small scales, M ~ 10~ 6 M Q . The abundance of the small- 
est subhalos can be inferred by extrapolating the subhalo 
mass function measured in simulations at larger mass 
scales. If we take the density within these 10 -6 Mq ob- 
jects to be ~ 100 times the smooth value, the radius of 
these objects is 10 -2 pc; the crossing time is 50 years; 
and the mean time between encounters is roughly a mil- 
lion years. For equilibration timescales less than the age 
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of the Sun (r© -C i©), the signal will be roughly at the 
equilibrium value of the smooth component for most of 
the time. However, for long equilibration timescales (e.g. 
the Earth) , the amount of depletion between interactions 
is negligible. This effect leads to a continuous build-up of 
WIMPs in the Earth, augmented by brief periods of an 
increased capture due to interactions with the subhalos. 
This results in an energetic-neutrino signal that today is 
higher than the signal that would be obtained from the 
smooth component. This can be understood as follows: 
For N <C N cq , the second term in Eq. ^ is small. While 
the cross section for the Solar System trajectory to inter- 
sect subhalos is cx /3 -2 / 3 , the capture rate while in them 
is oc /3, thus giving rise to a net increase in the capture 
rate oc /3 1 / 3 . Fig.[4]shows the net effect of this speed-up. 

In summary, we considered the effects of Galactic sub- 
structure on energetic neutrinos from WIMP annihila- 
tion in the Sun and the Earth. While DD experiments 
depend on the local dark-matter density and velocity dis- 
tribution, the energetic-neutrino fluxes from the Sun and 
the Earth depend on the past trajectory of the Solar Sys- 
tem through the clumpy Galactic halo. If experimental 
DD and vlT) signals are obtained before the dark-matter 
particle-physics parameters are known, then the poten- 
tial for probing dark matter via the DD / i/ID ratios will 
be compromised by the particle-physics uncertainties. If, 
however, the particle-physics parameters are known, then 
measurement of the DD rate and the vTD rates from the 
Sun and Earth can be used to probe Galactic substruc- 
ture at the solar radius, with the equilibration timescales 
setting roughly the mass scales that can be probed. 
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